Energy concentration and brittle crack 
propagation 



Abstract 

The purpose of this paper is to fill the gap between the classical treat- 
ment of brittle fracture mechanics and the new idea of considering the 
crack evolution as a free discontinuity problem. Griffith and Irwin crite- 
rions of crack propagation are studied and transformed in order to be no 
longer dependent on any prescription of the geometry of the crack during 
its evolution. The inequality contained in theorem 6.1. represents the link 
between generalized Irwin and Griffith criterions of brittle crack propa- 
gation. The physical meaning of this inequality is explained in the last 
section. 
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Note to the reader 



This is a paper which has been submitted for publication in a well-known jour- 
nal in 1997. The referee did not want to accept the paper unless substantial 
modifications are made. In my opinion the suggested modifications (except- 
ing ortographic, typographical or minor mathematical ones) were against the 
philosophy of the paper, which is: brittle fracture propagation is a geometrical 
evolution problem, therefore geometrical treatment is highly significant, both 
from mechanical and mathematical point of view. 

As this (eventually unpublished) paper circulated in manuscript, I think it 
is safe to depose it in the arXiv, for further reference and comparison with other 
works in this narrow domain. 

No modifications were made to the original file, excepting some preprint 
references in the bibliography, which were published since. 

1 Introduction 

Brittle fracture mechanics studies the evolution of the cracks in elastic bod- 
ies. Since the medium under consideration is supposed to remain elastic, the 
problem of brittle crack propagation concerns the evolution of the geometric 
support of the crack. We are facing here a problem placed in between geometry 
and mechanics. 

The field of research in brittle fracture mechanics has appeared in 1920 with 
Griffith's paper [G]. Among the basic references in this field we find: Irwin [I], 
Eshelby [Es], Gurtin [Gul], [Gu2]. 

Typically for classical fracture mechanics is that the geometrical nature of 
brittle crack evolution problem is obscured by the assumption of prescribed 
geometry of the crack. There are very few papers which do not make this 
assumption; as an example we cite here Ohtsuka [Ohtl — 4] and Stumpf & Le 
[StLe]. 

In the last decade a new current of ideas has emerged, starting with the 
article of Mumford & Shah [MS], submitted for publication in 1986. It is the 
first time, to our knowledge, when the crack itself is considered as the unknown 
of the problem of brittle crack appearance. Some new mathematical results of 
De Giorgi and Ambrosio ([DGA], [Al], [A2]) in the field of geometric measure 
theory, set the energetic model of brittle crack appearance proposed by Mum- 
ford and Shah in the functional space SBV, of special functions with bounded 
variation. 

The strange unknown of the problem — the crack — is replaced by a more 
familiar one: the displacement field of the cracked body. This field is allowed 
to be discontinuous, being a special function with bounded variation. Recent 
papers study the space SBD of special functions with bounded deformation, 
as Bcllcttini, Coscia & Dal Maso [BCDM] and Ambrosio, Coscia & Dal Maso 
[ACDM] , naturally associated with the expression of free energy potential of an 
elastic body suffering small deformations. 

It is now clear that a promising point of view concerning crack evolution 
is to consider it as a free discontinuity problem. The purpose of this paper is 



2 



to fill the gap between the classical treatment of brittle fracture mechanics and 
these new ideas. We present here some of the results obtained in [Bu]. The 
Griffith and Irwin criterions of crack propagation are studied and transformed 
in order to be no longer dependent on any prescription of the geometry of the 
crack during its evolution. 

The content of the paper is described further. In the section "Preliminaries 
and notations" some basic constitutive assumptions are made. The following 
section, " Statics of a fractured body" , contains a mathematical treatment of the 
equilibrium of an elastic cracked body in a functional setting compatible with 
the space SBV. This section contains also a brief description of the Dirichlet- 
to-Neumann map of a cracked body. For more information about this notion 
we send the reader to [C], [SU1,2]. 

The section " Kinematics of crack propagation" starts with the introduction 
of a set of reasonably smooth endomorphisms of f2, the reference configuration 
of the body useful in the sequel. The notion of smooth crack propagation is 
introduced, in various forms (definitions 4.1. to 4.3.). A smooth crack propaga- 
tion can be seen as a curve of endomorphisms of Q, namely 1 1— ► <f> t . The initial 
crack K, which is a surface with boundary in O, becomes at the moment t the 
actual crack K t — <j> t (K). We mention that in [StLe] variations of the crack by 
diffeomorphisms arc considered; also, in [Ohtl — 4], a curve of diffeomorphisms 
1 1 ► is associated to a curve of cracks 1 1— ► K t . 

In the section " Crack propagation criterions" Griffith and Irwin criterions of 
brittle crack propagation are reformulated for smooth crack propagation. This is 
accomplished by proposition 5.1., definitions 5.1. (generalized Rice's J integral) 
and 5.2 (concentration coefficients of the clastic energy). 

The section "Energy concentration and crack propagation" begins with the 
introduction of the measure | K2 | (u) (definition 6.1.) associated to a field u 
of displacements. This measure describes the distribution in the body of the 
energy release rate due to the crack propagation. In definition 6.2. is introduced, 
for the same arbitrary field of displacements u, the upper energy concentration 
coefficient as a measure, namely CM + (u). The main result of the paper is 
theorem 6.1., which describes the relations between these two measures. 

A byproduct of the proof of theorem 6.1., with important physical signifi- 
cance, is presented in the last section. 

2 Preliminaries and notations 

Let ft C R n represent the reference configuration of an elastic body. The 
space dimension n equals 2 or 3 and fl is an open bounded set with piece- 
wise smooth boundary dfl. The volume density of the free energy is a function 
w = w(Vu). u denotes the displacement of the body from the reference config- 
uration. We suppose that the function w depends only on e(u), the symmetric 
part of Vu. In this paper we work with a quadratic expression of w: 

io(Vu) = ^CVu : Vu . (1) 
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The tensor C has the symmetries: 

Cijkl = Cjikl — Cklij ■ (2) 

Most of the facts herein are true whenever w is a convex, C 2 function, satisfying 
the following growth condition: there arc two positive constants c, C such that 

\/E 6 R nxn , c | E sym | 2 < w{E) <C\E\ 2 (^E sym := ^{E + E T ) 

' (3) 

We will restrict however our attention to the case when the energy function w 
satisfies ©, has the form (JTJ and C satisfies (J2J. 
The first Piola-Kirchoff stress tensor is defined by: 

dm 

<?ij = o-ij(Vu) = ^tt(Vu) 

The body evolves in a quasi-static manner, in the absence of volume forces, so 
at any moment the stress tensor a is divergence free. 

The Hausdorff fc-dimensional measure is denoted by U k . H 2 is the area 
measure and Ti. 1 is the length measure. In order to simplify the denominations 
we shall often call TU 1 ^ 1 the area measure. The volume measure is H n = 
the Lebesgue measure. 

Any crack in the body is seen as a "crack set" (see Ball [Ba]): a closed 
countably rectifiable set K C f2 with finite area. Any smooth hyper-surface 
with bounded mean curvature and with boundary in Q is called a "smooth 
crack set". A smooth crack set it is therefore a crack set which carries the 
geometric structure of a manifold with boundary. The boundary of the smooth 
crack set K, denoted by dK, represents the edge of the smooth crack set. The 
edge dK is a n — 2 piecewise smooth surface without boundary. 

In this paper we work with smooth crack sets K with the property: 

(Jf n 0fl) = o . 

This assumption has the following meaning: the intersection of the smooth crack 
set with the exterior boundary of the body dQ is at most a reunion of curves. 

The space SBV(f2, R n ) of special functions with bounded variation was in- 
troduced by De Giorgi and Ambrosio in the study of a class of free discontinuity 
problems ([DGA], [Al], [A2]). For any function u e L l {£l,R n ) let us denote by 
Du the distributional derivative of u seen as a vector measure. The variation 
of Du is a scalar measure defined like this: for any Borel measurable subset B 
of the variation of Du over B is 

| Du | (B) = sup |^ | Du{Ai) | : At C B , A, n A j = Vz ^ j j . 

The Lebesgue set of u is the set of points where u has approximate limit. The 
complementary set is a C n negligible set denoted by Su- 
The space SBV(f2, R n ) is defined as follows: 

SBV(0, R n ) = {u G L\n, R n ) : \ Du \ (O) < +oo , | D s u | (fi \ S u ) - 0} . 
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Let us define the following Sobolev space associated to the crack set K (see 
[ABF]): 

W 1 / = ju e SBV(fi,iZ") : J I Vu | 2 dx + J [u] 2 dU^ 1 < +00 , | D s u |< H\ 
It has been proved in [DGCL] the following equality: 

w^ 2 (n, R n ) n J?™) = w 1 ' 2 ^ \ k, R n ) n L°°(n, R n ) . 

We shall need further a norm with physical dimension on the Sobolev space 

W p ' 2 '°°(Cl, R n ) = W p < 2 (fl, R n ) n L°°(Cl, R n ) . 

Consider therefore A as an unit of length. We take the following norm on 
W p > 2 (n,R n ) : 

||u|| p , 2 = ||u|| L2 + A ||Vu|| i2 + ... + X p \\V p u\\ L 2 . 
This norm induces on W P ' 2 '°°(Q, R n ) the norm 

||u|| p ,2,oo = max(vol(Cl) ||u||oo, ||u|| p , 2 ) , 
where vol(Cl) is the volume of SI measured with the unit A". 

3 Statics of a fractured body 

3.1 Functional assumptions 

Let us consider an elastic body Q with a smooth crack set K . The displace- 
ment u° G H2(dQ,R l ) fl L°°(dQ, R n ) is imposed on the exterior boundary 

an. 

Therefore the equilibrium displacement u is a special function with bounded 
variation. ^.Frorn the Calderon & Zygmund [CZ] decomposition theorem we 
obtain the following expression of Du, the distributional derivative of u seen as 
a measure: 

Du = Vu(x) dx + [u] O n dH^ 1 . 
We deduce from here the Stokes formula: 

{Du,<t>) = f Vu-4>dx+ [ [u] ® n • dH^ 1 V0 e Cg°(fi, R nxn ) . (4) 

This formula is the mathematical expression of the fact that are no concentrated 
forces on the edge of the crack. 

Any element of W^ 2 , compatible with the boundary condition, is called an 
admissible displacement. The family of admissible stresses will be defined in 
the sequel. 
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The jump of a field a across K is [a] = a + — a . Del Piero considers in [DP] 
the following set of admissible stresses for a fractured media: 

Y(fl) = {a = cr T e L 2 (n,R nxn ) : diva = 0} . 

He proves that if a G Y(Q) then er+n = a~n. This fact allows us to use the 
notation an without confusion. 

We connect further the definition of an admissible displacement field with 
the definition of an admissible stress field. 

Definition 3.1.1. Let u° e (dQ, R n )nL°°(dn, R n ) represent an imposed 
displacement on the exterior boundary of the body fl and let K be a crack set in 
il. The set of admissible displacements with respect to u° and K is 

W£ 2 '°°(u ) = {u eW^ 2 (n\ K, R n ) n L°°(n, R n ) : u = u° on dfl } (5) 

and the set of admissible stresses with respect to K is 

Y K (fl) = {a = a T G L 2 (n, R nxn ) : diva = , an = on K] . (6) 

On the crack set the body is force free, i.e. <r + '~n = on K 1 where n is the 
field of normals on K and a + , a~ are the traces of a on the sides of K. The 
sign convention is taken such that the field of normals to K, n, points to the 
"+" side of K. There are no concentrated forces on the edge of the crack. 

3.2 Dirichlet-to Neumann map 

The equilibrium displacement of the body, when the boundary displacement 
u° is imposed, is solution of the problem: 

div cr(Vu) =0 in Q \ K 

cr(Vu)+n = cr(Vu)-n = on K (7) 
u = u° on <9fi 

The equilibrium displacement u minimizes the functional 

E{\) = / io(Vv) dx 
Jn 

defined over all v G W^ 2 (Q \ K, R n ), v = u° on dtt. The solution of the 
minimization problem is unique to an arbitrary piecewise rigid displacement 
equal to on dil. Also, u is smooth and essentially bounded. 

Yk{&) is a subset of Y(Q); in the definition the quantity an does not depend 
on the choice of n. 

The polar of the convex functional W{e) = J n w(e) dx is (see Moreau [M]): 
W*(a) = sup {(a, 6} - W(e) : e G L 2 {^ R^)} . 
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For any admissible stress a € Yk (O) and for any minimizer u of the func- 
tional E over the class of admissible displacements, the following inequality is 
true: 

f (an) • u dH n - 1 - W*{a) < E{u) . (8) 
Jdn 

The inequality becomes an equality if and only if a = er(Vu). 



For a given crack set K the Dirichlet-to-Neumann (or the response) map 
associated to the domain ft with the crack K inside can be defined. This map 
is defined in the following way: let u° be given on the boundary of f2. The 
solution of the Dirichlet problem 

div cr(Vu) =0 in D, \ K 

u = u° ondVt (9) 

er(Vu)n = on (both sides of) K 

will be denoted by u(K, u°). This solution is uniquely determined to a rigid 
displacement equal to on dfl hence <r(K, u°) is unique. 

It is obvious that u(K, u°) is also solution of the Neumann problem 

div cr(Vu) = in Q \ K 

cr(Vu)n = <r{Vu(K,u°))n on dQ (10) 
er(Vu)n = on (both sides of) K . 

Therefore the function u° i— > <r(Vu(K, u°))-n maps naturally a Dirichlet bound- 
ary condition to a Neumann boundary condition. The Dirichlet-to-Neumann 
map is the linear application: 

T(K) : H^(dn) -> H-i(dSl) , 

(T(if)u°,v°) = f a(Vu(K,u°))n-v dH 11 - 1 . (11) 

Here ( , ) means the duality product. The Dirichlet-to-Neumann map is con- 
tinuous due to the well known continuous dependence of the solutions of the 
problems © and (|10fl with respect to boundary data. The symmetries of the 
elasticity tensor C make this map self adjoint. 

If u° is essentially bounded then the solution of the problem (unique to a 
rigid displacement equal to on 90) 

minimize J fi w(Vu) dx over all u € W^ 2 ' 00 ^ ) 
is u(K , u°), hence 

i (T(^)u°,u ) = min w(Vu) dx : u6^ 2 '°°(u ) 
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4 Kinematics of crack propagation 



Let us consider a hyper-elastic body with the reference configuration f2. At 
the moment t — in the body there is the crack set K. The evolution of the 
crack is a curve t i— » K t such that Kq — K and for any t < t' K t C Kf ■ 

We are searching for an easy way to describe the evolution of the crack 
set. If the dimension of the space is n = 2 then an obvious way to describe 
the evolution of the crack set is to consider only the evolution of the edge of 
if, i.e. dK, which is formed by a finite number of points. In the general case 
n = 3 this corresponds to the evolution 1 1— > dK t where dK t is a curve. Because 
the evolution of a surface in the three-dimensional space is conceptually simpler 
than the evolution of a curve, the before mentioned choice of modeling the crack 
propagation does not simplify our problem. We prefer instead to think at the 
crack set K t as a deformation of the initial crack set if by a diffeomorphism 
(f>t : ft — ► Cl. We replace therefore the curve t i— » K t with the curve t i— > (fit, 
satisfying the condition 4>t{K) = K t . (However, not any curve t i— > K t has an 
associated curve t i— ► (fit- As a counterexample we may think at cracks which 
after a time develop new branches.) 

The curve 1 1— ► <f>t lies in a set of diffeomorphisms. 

Let us consider the following family of diffeomorphisms: 

M c = {(/.eC^^nW^' 2 ' 00 ^,^ 1 ) : (12) 

is a diffeomorphism and supp (</> — 1q) C } 

The set M c can be conveniently completed to a topological group, which has 
a geometrical structure. 

We consider, inspired from Ebin & Marsden [EbM], the closure of the set 
M c in the W s ' 2 '°° topology: 

D s /-°°(n) = Cl ws ' 2 '°° M c , (13) 

where s is chosen to be greater than a critical value. 

The tangent space at D s ' 2 >°°(n) in l n is denoted by T ln D s ' 2 ^(Q) and it is 
equal to 

Wo' 2 '°°(n,R n ) = Cl WS ' 2 '°° {ry e C°°(Q,R n ) r\W s ^-°°{n,R n ) : supp r? C fi} . 

This tangent space may be endowed with the norm || ||i s2 = | ||vj/s,2,«.. 
The tangent space at <p and the corresponding norm can be obtained from the 
tangent space and norm at the identity by right multiplication with <j>. 

The existence of the exponential map is proven in [EbM], theorem 3.1. for the 
group of W s ' 2 diffeomorphisms of M, where M is a compact manifold without 
boundary. The group defined at l|13H can be seen as a group of diffeomorphisms 
of a compact manifold equal to identity near a point of the manifold. The 
results of Ebin & Marsden are still true in this case. In theorem 3.4. from the 
same article is proven that any continuous time dependent vector field can be 
integrated. We give the exact statements of these theorems, adapted to our 
case. We do not give the proof, since it requires the same techniques as the ones 
employed in the article mentioned before: 
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Theorem 4.1. Let Q, be an open bounded set in R n with piecewise smooth 
boundary and s a number satisfying 



Let i] G Wq 2,c °(VL, R n ) be a vector field on 0, vanishing near dVt. Then the 
one parameter flow generated by rj is a C 1 one parameter subgroup of D S C > 2,0 °{£1). 

Let 1 1 — ► Tjt G Wo' 2 (f2, R n ) be a continuous time dependent vector field on il. 
Then the problem: 

Vt = 4>t4^ X , 0o = In 
has a solution, unique, which is a C 1 curve in _D^ 2, °°(r2). 



We shall define further the notion of smooth crack propagation curve. 

Definition 4.1. Let us consider s > f + 2. An s-smooth crack propagation 
curve is a C 1 curve 

te[0,T]^&eD s /' co (n) , 

satisfying the following properties: 

• 0o = In, the identity map ofVL, 

• 4> t {K) C 4>t'(K) for anyO<t<t'<T. 

For any smooth crack propagation curve we associate the curve t K t = 
(j>t(K). The last curve has the obvious property that for any t > the crack 
set K t can be continuously deformed into K. In other words, K t and K are 
topologically the same. We restrict our attention to this kind of evolution of 
the initial crack set. In this approach the initial crack K may be as complex as 
we wish, but this complexity remains the same during its evolution. 

There are infinitely many smooth crack propagation curves 1 1— ► <j>t with the 
same associated function 1 1— » K t . 

Definition 4.2. Let t e [0,T] ^ u°(t) e C(dQ,R n ) be a C 1 curve of 
imposed displacements on the exterior boundary of the body. An admissible 
fracture curve is any C 1 function 

t € [0,T] i ► «> t ) G W 1 / x D 2 C P-°°(Q) , 

satisfying the following items: 

• for any t G [0, T] u* = u°(t) on dfl, 

• t i— > <j> t is a s-smooth crack propagation curve. 
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Let us denote u t — ujf.^ . Because for any t <j) t equals the identity map 
near the boundary of 0, we deduce that u t = u°(£) on dfl. Remark also that 
the application u* i— > u*.^ -1 , where <f> G £)* ,2: °°(f7) , maps W^ 2 onto W^ K y 
therefore u t is an admissible displacement of the body with the crack set 
UK). 

Definition 4.3. A balanced fracture curve is an admissible fracture curve 
t i ► (uj,0t) suc h that at any moment t>0 the displacement u t = u£.^~ is a 
solution of the problem 

mm | jT »(Vu) d,:ue W$"(u°(t)) J . 



For any smooth crack propagation curve there is an associated balanced 
fracture curve, unique to rigid displacements. 

We are concerned now with the evolution of the area of the crack set. 

It is well-known (see [All]) that the variation of the area of the crack set 
UK) is 



±-H n -\MK)) = I divsrn dTT- 1 

dt Ht(K) 



where div s is the tangential divergence with respect to the surface 4>t (K) . If we 
denote by n the field of normals to this surface, the expression of the tangential 
derivative div s is: 

div s ri = ~ //,., 11,11, . 

For a smooth crack propagation curve t <j) t the condition that the crack 
grows implies that for any t > 

Vt-n =0 on UK) , 

S MK) div sm <m n ~ l >o w 

It will be useful to define a perimeter measure of the edge of the crack. For 
this let us consider a crack set K. The perimeter of dK can be defined, with 
the help of a flux-divergence formula, like this: 

T(dK) = sup^J div s r] ATT^ 1 ; i\ G C °°(O,i?") , | rj |< 1 , j) • n = on if J . 

The perimeter of dK as a measure is defined first as an additive function over 
the set of open subsets B of Q: 

T(dK)(B) = sup^J div sV dTT" 1 ; rj G C^{B,R n ) ,| rj |< 1 , rj • n = on 

(15) 

This function generates a measure which is the perimeter of dK as a measure. 
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5 Crack propagation criterions 



In order to select one or more crack propagation curves among all admissible 
ones one needs a criterion of brittle crack propagation. There are several such 
criterions. We shall discuss about two of them, in principle different. 

The Griffith criterion asserts that during the crack propagation the energy 
release rate due only to the crack evolution has to be greater than a critical 
quantity. This criterion may take different mathematical forms. Whatever this 
form may be, this criterion is formulated in terms of a critical speed. 

The Irwin type criterion asserts that during the crack propagation some 
intensity (or concentration) factors have to be greater than a critical value. 
The original Irwin criterion is formulated in terms of stress intensity factors. 
However, it is straightforward that it can be reformulated by using the elastic 
energy concentration on the edge of the crack. We shall refer to any crack 
propagation criterion using energy concentration factors (whatever they mean) 
as to the Irwin criterion. 

The goal of this section is to give precise mathematical expressions to the 
Griffith and Irwin criterions for the class of smooth crack propagation curves. 

We shall use in the sequel the notation 

u(0t,u°(t)) = u t = u(MK),u°(t)) 

since the initial crack K is given. The stress field associated to this displacement 
is a t = cr(Vu t ). 

The power communicated by the universe to the body at the moment t has 
the expression 

P(t) = [ a t n- u° dH' 1 - 1 . 
The Griffith criterion of brittle fracture propagation is the following: 



A balanced fracture curve t i— > (u*,<^ t ) satisfies the Griffith 
criterion if at any moment t > 



^|^(Vu(^,u°(i))) dx + GK' 



■ t (K))\<P(t) . (16) 



G is a material constant, named the constant of Griffith. 

The energy release rate due only to the propagation of the crack has the 
well-known expression: 

£(t) = P(t) - j t jf«;(Vu(0 t) u o (i)))cLc . (17) 

The inequality i|16|) is equivalent to 

£(t) > g i. tt-HMk)) ■ 
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This criterion seems hard to work with it. The following proposition will 
lead us to an easier version of l|16|) . 



Proposition 5.1. Let t i— » fa be a smooth crack propagation curve, r\ t = 
fa-fa^ 1 and u° G H*(d£l) nL°°(dfl). Then for allt>0 the following inequality 
is true : 

(|[T(0 t )]u°,u°) < (18) 



< 



{[CVu(0 t ,u°) : Vu(0 t ,u )] div tft - 2 [CVu(^, u )] .. [Vu(&,u°)]. fc [V % ] fej } dx 



Proof: The smooth crack propagation curve fa is associated to the balanced 
fracture curve (u(fa, u°).^ t , fa). Let us choose a t > and keep it fixed. We 
define now, for any s > 0, 

u ts = u(fa,u°).fa.faj+ s . 
It is obvious that for all s > S Ut s C fa +s (K), hence 

Vs>0, (T(0 t+s )u°,u o ) < f CVu t , s : Vu u dx . 

For s = u tl0 = u(fa, u°), so 

(T(^)u°,u ) - / CVu f , : Vu t , dx . 



We deduce that 
(- 



[T(0 t )]u o ,u o } = -^-(T(0 t+s )u o ,u o ) u _ o < A / CVu, s : Vu, s dx u _ . 
ds 1 ds J n 

It remains to prove that 

{[CVu(&,u°) : Vu(0 t ,u°)] div T) t - 2 [CVu(&, u )] .. [Vu(&,u°)] <fc [V^} dx 
d /" 

= — J CVu tiS : Vu M dx| s=0 . 

The last equality results by direct calculation. It is sufficient to make in the 
integral 

CViif s : Vu f , s dx 



the change of variables y — fa.fa, s [x) and then to perform the calculation of 
the derivative. • 

Definition 5.1. The generalized Rice's J integral is the functional 
K2 : D s /'°°{fl) x T ln D s c ' 2 '°°(n) x nL 00 ^)} -> R , 
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K2(4 > , V ,u ) = j- 1 [CVu(0,u°) : Vu(0,u )] div 77 + (19) 
+ [CVu(^u o )] y .[Vu(0,u°)] tfc [Vr,] fcj } dx . 

We shall explain why K2 is called "the generalized Rice's J integral" in the 
remark 5.2.. Before that we give an alternative expression of K2. 
Let us simplify the notations: 

o- = CVuO,u ) , u = u(0,u°) , w = icVuO,u°) : Vu(0,u°) . 

If u is a C 2 function then the following equality makes sense: 

W T)i,i-VijUi,k Vk,j = [w T)i ~ 0"yUj ifc IJ^j-fftaU^i Vi+^li^lM Vk+^li,i^l,k Vk ■ 

The divergence of the field <r is equal to 0, hence 

#2(0, 77,11°) = - [w m - (T U n Uk rj k ] . dx . 
Jn 

Let d(j){K) be the edge of the crack set <p{K). We define the tubular neighbor- 
hood of d(j)(K) of radius r: 

B r = B r (d<f>(K)) = U xedHK) B(x, r) . 

The field of normals over d B r (d(f>(K)) will be denoted by v, without specifying 
the parameter r. 
The function 

x e fl w divr](x) - a H u lyk rj k ,i (x) 
is measurable, 77 = on <9f2 and an = on <fi(K), therefore 

K2{<j>, 77, u°) = lim \ f {-tori ■ v + <T U Mi, k Vk »i\ dH^ 1 U / [10] 77-11 dU 

[J dB r (d<P(K)) J J<t>(K) 

If ?7 • n = on <j>(K) then 
K 2(0, 77, u°) = lim J / - {^77 • ^ - oum, k ri k Vi} dU^ 1 \ . (20) 

[J dB r (d<j>{K)) J 



Remark 5.1. In [Ohtl — 4] Ohtsuka generalizes the Rice's integral J. For a 
given, very smooth, curve of increasing crack sets t ^ K tl he founds a vector 
field 770 whose one-parameter flow t has the property 

K t = MKo) ■ 



n—1 
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He takes into consideration not only imposed displacements on the exterior 
boundary of the body but also imposed forces. The connection between his 
generalization and the functional K2 is the following: 



GJ(K t ,u°) = K2{0 t , m ,u°) . 
Ohtsuka proves that GJ(Kt,u°(t)) equals the energy release rate £(t). • 

Remark 5.2. In the very particular case of one-dimensional crack prop- 
agation following a straight line Rice [R] expresses the energy release rate by 
an integral, named J. The expression of J is formally similar with the right 
member of (12(H when 77 is taken equal to a constant vector (the speed of the 
crack edge) parallel with the crack. 

The source of the notation K2 can be found in the particular case of a 
straight crack in the 2-dimensional configuration. Let us suppose that the body 
is a cylinder with section Q and consider only anti-plane displacements: 

u : Q x R — > R 3 , u(xi,X2,xs) = (0, 0, u(x%, X2)) ■ 

In this case the integral J is proportional with the square of the stress intensity 
factor in mode III denoted by Km. • 

In the sequel we shall use the notation: 

Tfo) = T(4>(K)) . 
The energy release rate £ (t) has the following expression: 

£{t) = (T&)u (t),u°(t)) - \ |{(T^)u°(t),u°(t)}} . 
We use the fact that T(</>) is self-adjoint for proving that 

<T0 t )u o (i),u°(i)} = (T<p t )u°(t),u°(t)) . 
The last two relations show that the energy release rate has the expression: 

m = -^|[T(^)]u W,u°W) • 
We use proposition 5.1. and definition 5.1. for proving that 

£{t) > K2(<j> t ,r) t ,u°(t)) . (21) 
We propose now a Griffith type criterion of smooth fracture propagation. 

A smooth crack propagation curve t 1— * <f> t satisfies the gener- 
alized Griffith criterion if at any moment t > 

K2((b u r lt ,u°(t)) > G \\rit\\oc'P(d4>t(K)){supp rjt) (22) 
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In our criterion of brittle fracture propagation £(t) is replaced by K2{<f>t, r]t, u°(t)) 
and the variation of the area of the crack 

/ div s r, t dH' 1 - 1 
JMK) 

is replaced by \\r] t \\oo'P(d(j>t(K))(supp if). This criterion is therefore stronger 
than the classical Griffith criterion (this is a consequence of l|21[l and of the defi- 
nition of the perimeter (|15fl ). The two selection criterions are however equivalent 
if at any moment t the crack is smooth enough and it evolves in a very smooth 
manner. 



The Irwin type criterion of brittle crack propagation is formulated in terms 
of concentration coefficients of the elastic energy. 

Definition 5.2. Let u be a special function with bounded variation. The 

lower (respectively upper) concentration coefficients of the elastic energy of u 
are 

C2~(u, x) = - lim inf JnnB ^ _ , (23) 

7T r^O r" 

C2+ (u, x) = - lim sup J " nB( *' r) — . (24) 

7T r —i0 T n 



If x € fl \ Su then C2 + (u, x) = 0, because u>(Vu) has approximate limit in 

x. 

Let us denote by Eu the symmetric part of the distributional derivative of 

u: 

Eu = - (Du + D T u) . 

Eu can be seen, like Du, as a vector measure; let | Eu \ be the scalar measure 
of the variation of Eu. For any borelian set B C il 

| Eu | (B) = sup jf] | EuiAi) | : u£ t A* C B , A t n Aj = V* ^ i| • 

Let us consider the set 6u introduced by Kohn [K]: 

f ^ , I iJu I (B(x,r)) 
Q u = <xeQ: limsup J 1 v _) - " > 

A result from [ACDM] , concerning special functions with bounded deformations 
in the particular case of u special function with bounded variation, assures us 
that 0u differs from Su by a set which has W n_1 measure zero. We use the 
growth condition of w 10 in order to prove that if x S 0u then C2 + (u, x) = 
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0. Therefore we have the following localization of the support of the energy 
concentration coefficients: 

suppC2~(u, ■) C supp C2+(u, •) C S u \ O u ■ 

One can prove, using the same proof, that if if is a crack set and u € W^ 2 fl 
L°°(n,R n ) then 

supp C2 _ (u, •) C supp C2+(u, •) C dK . 
The Irwin type criterion that we propose is the following: 

A smooth crack propagation curve t i— > <j> t satisfies the Irwin 
type criterion if at any moment t >0 and for any x G dipt (K) 

C2+(u(0 tj u°(t)),aO > G . (25) 



6 Energy concentration and crack propagation 

The goal of this section is to find the connections between the Griffith cri- 
terion (1221) and the Irwin criterion (|25f) . using a measure-theoretic approach. 
The first step is to associate a measure to the generalized Rice integral if 2. 

Definition 6.1. Let us consider the linear functional 

k2 : SBV(fi,iT) ni°°(fi,i? n ) x C* °°(rj,i?") , 

fc2(u,7?) = / {o-(Vu)ijU itk rik,j - ffl(Vu) div 77} cfc . (26) 
Ju 

Then, for any u e SBV(^, i?™) n L°°(fi, R n ) and for any open subset B of n 
we define 

I if2 I (u)(B) = sup{k2(u, rj) : supp 77 CC B , | 77 |< 1 , rj ■ n = on Su} • 

(27) 

Since \ if 2 | (u)(-) is finitely additive, it gives raise to a positive measure, 
denoted by \ if 2 | (u) too, named the generalized Rice integral as a measure. 

It is easy to see that for a smooth crack propagation curve 1 1— ► (f>t it is true 
that 

Vft>0, if2(^,7j t ,u (t)) = k2(u(4> t ,u°(t)),r h ) < \ K2 \ (u(^, u°(t)))(supp Vt ) 

Remark 6.1. The measures V(K) and | if2 | (u) are constructed in the 
same way. The departure point was in both cases a linear functional over 
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Cq°(Q, R n ) — the area variation integral in the first case and the k2 integral in 
the second case. • 

We shall define now the upper energy concentration coefficient as a measure. 

Definition 6.2. Let us consider u e SBV(f2, R n ). For any r > let B r (u) 
be a tubular neighborhood o/Su of radius r. For any open set B in Q we define 
the upper energy concentration coefficient: 

CM+(u)(B) = limsup JB ^ nB . (28) 

Since CM + (u)(-) is an additive function, it give raise to a positive measure 
denoted by CM + (u) too, named the upper energy concentration coefficient as a 
measure. 

CM + (u) splits in two parts: the absolute continuous part with respect to 
C n and the singular one. 

CM+(u) = CM+(u) A + CM+(u) s . 

Let us consider, for u e SBV(£1, R n ), the decomposition of flu: 

Du = vudi + [u] ®n an?- 1 . 

The measure CM + (u) s can be decomposed in two parts — the absolute con- 
tinuous part with respect to the measure 



J(u) = [u](8)ndWp 



jn-l 

'u 

and the remaining singular part — i.e. : 

CM + (u) s = CM + (u) J + CM + (u) c . 

It is easy to see that CM + (u) A = 0. We obtain therefore the following 
decomposition of CM + (u): 

CM+(u) = CM+(u) J + CM+(u) c (29) 

Ifue W^^Su.^and S u is a crack set, then we find, by direct calculation, 
the expression of CM + (u) 



c. 



CM+(uf = C2+(u,x) &H\ 1 - 2 . (30) 

aS u 



The main result of the paper is the following: 
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Theorem 6.1. Let K be a crack set and u 6 w£ 2 '°°(u°) R w£' c 2 (ft\ AT, R n ). 
Then u = u(A", u ) if and only if for any open set ficSl 

| A2 | (u)(B) < CM+(uf(B) . (31) 



A consequence of theorem 6.1. and (|31[1 is that the Griffith criterion (1221 is 
stronger than the Irwin type criterion i|25|) . Indeed, let us consider a smooth 
crack propagation curve t i— » (fit satisfying the generalized Griffith criterion (|22|l . 
Then, at any moment t > 0, we have: 

| K2 | (u(&,u°(f)))(fi) Halloo > A2(</) t ,77 t ,u (t)) > G||r ?t || oo P(a0 t (A'))( SU ppr ?t ) . 
This chain of inequalities, l|31|) and (|3U[1 imply the following relation: 

/ C2 + (u(0 t ,u°(i)),.) dTT- 1 > G V(dMK))(supp rn) . 

JdMK) 

Therefore at any moment t > there is at least an x € d(fit{K) such that 

C2 + {u((fi t ,u°(t)),x) >G . 



7 Proof of theorem 6.1. 

1 st implication. We want to prove that if u = u(A, u°) then for any open 
set B the relation (|31|) is true, namely 

| K2 | (u)(B) < CAf+(u) c (B) . 

Let us consider a vector field 77, such that: 

r) G C °°(f7,A") 

' " 'n a ( 32 ) 

77 • n = on A 

supp n d B 

Any such vector field 77 generates an one-parameter flow s >— > <fi s . We can always 
find a curvilinear coordinate system {ot\, ...,a n -x,j) such that on the edge of 
the crack set dK 7 = 0, the surface 7 = s is the boundary of a tubular 
neighborhood of d A, namely 

B r ( s )(dK) = U xedK B r{s} (x) 

and for all s > 

s (A)c KUB r(s) (dK) . (33) 
We can suppose moreover that 

lim^l = 1 

s->0 s 
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because (|32H affirms that the crack 4> S (K) evoives with sub- unitary speed. 
We denote by u s the soiution of the problem: 

div cr(Vu) = in O \ (K U -B r (s)) 
a(Vu)n = on dB r{s) (dK) " 

cr(Vu)n = on (both sides of) K\B r[s) (dK) [ ' 

u = u° on <9f2 . 

We consider, for any s > 0, the following stress held: 

*.W = I n (VUs)W ^\ n \( ifU ^w) (35) 
w [0 elsewhere . v ; 

It is obvious that for any s a s is admissible with respect to <f> s (K) (definition 
3.1.1). We use © and proposition 5.1. to deduce that 

L i«(Vu s ) dx — f n WVu) dx , , 

liminf^— i ^— 1 < 36 

s^O s 

L w(X7u(d) s , u )) dx — f n w(Vu) di: , n% 

< lim^— i — — ^— < -A'2 1 o ,?7,u . 

s^Q s 

Our goal is to prove that 

,. . c Lw{Vu s ) dx - Lw(Vu) dx ,. L,nii m ( Vu ) dx 
hm mf ^ — ^ — = - hm sup — — — . 

s->o s s ^o r(s) 

(37) 

For this it is sufficient to show that 

fo\ a w(Vu s ) da; — L, „ hi(Vu) da; 
hm M ^> MBr ^ = . (38) 

s^O s 

Indeed, if l|38() is true then an easy calculation shows that (|37|l is true. With 
this in mind we return to i|36|) and we see that 

(' nR io(Vu) da; 

_ limsup ^^' 05 . . < -K2(l a ,V,n°) , 

s^o r(s) 

hence 

f R nR w(Vu) da; 

We take the supremum of the left side term of the last inequality, over all ij 
satisfying (|32p. and we obtain the inequality 

I' o w(Vu) da; 

|A2|(u)(B) < limsup r(a) . . . 

s ^o r(s) 

This inequality is equivalent with (|31|l (as it is shown by definition 5.2. and i|30|) 
)• 
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Proof of i|38|) . Wc consider for any s the function 

M _ j u s {x) if x<£B r{s) (dK) 
s( > ~ \ Q it x e B r(s) (dK) . 

Each u s is a special function with bounded variation on 17. 

In the following we shall proceed as in [BBC]. We denote by ej., e n _i, e 7 
the local basis associated to the previous mentioned system of coordinates. 
e 1 , e n-1 ,e 7 represents the dual local basis. W can choose the system of 
coordinates such that e 7 = /in, where n is the normal of the surface 7 = ct. 
and h a scalar function. 

The expression of the divergence of a tensor field in this system of coordinates 

is 

div a = (V e „. a) ■ e Qj + (V e , a) ■ e 7 . 
We make the notation 

div T a = (Ve Q . cr) ■ e a% 

and we remark that in div T a enters only partial derivatives with respect to a,. 
The divergence of the field a can be written as 

Id. , 1 dn 
div a — div T a + —7— (cr • nj — — a ■ — . 

h oj h 07 

Consider a — cr(Vu) and 

a m (ai,s) = cr s (a,,s) . 

The field o~ m has the following interpretation: imagine that at the "moment" s 
the curvilinear cylinder B r ^(dK) is removed from the elastic body f2. Then 
cr m (-,s) represents the superficial stress of the remaining body on the new- 
created surface, i.e. dB r r s \(dK). Therefore 

a m (ai,a)n = , 

where n is the normal of the surface 7 = a. 

We consider, for any s the auxiliary elastic problem: 

div a^ ux = in Q \ (K U B r(s) {dK)) 

a aux n _ o-^sjn on (K\B r{s) (dK))udB r(s) (dK) (39) 
u a s ux = on 9fi 

with the solution u® ux (ai,j), 7 < s. The associated stress field is a^ ux (ai, 7). 
The superposition principle applied to the elastic problem ^ with the solution 
u = u(K, u°), H34|l and (|39|l leads us to the equality 

a( ai ,s) = a m (a l ,s)+a a s ux (a u s) . (40) 

The following equality is obtained from a linear combination of the divergence 
equations contained in the problems JHJ, (|34() and (|39|l , written in curvilinear 
coordinates, together with JJQJ. 

Q Q 

— (an) (at,7 = s) = — {a% ux n) (^,7 = s) - Hid<t>„7 = s) . (41) 
07 07 
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We make the derivative with respect to s in the equality (|4(J|) , then we multiply 
by n and we obtain 

^ (an) (a u 7 = s) = A (cx^n) 7 = «) + ^ (a? M n) (a 4 , 7 = *) . 

(42) 

(a7 x n) (a*, 7 = s) = -h div a m (ai,-y = s) . (43) 



lHH and g5J lead us to the equality 



We keep now s constant. Let us consider s' € [0, s\. We use the notation 
n s > = fi\ (K U B r(sl) (dK)) , 

S s , =dB r{s , ) {dK) . 

It is obvious that if s' < s" then f2 s » C f2 s / . 

The field a^ x is by definition the elastic stress field in the body tt s > , with the 
imposed displacement on it's exterior boundary dtt equal to and subjected to 
the loads a(cti, r(s'))n(ai, s') on it's interior boundary S s > ■ It is straightforward 
that the same field represents the elastic stress in the body il s , with the imposed 
displacement on it's exterior boundary dfl equal to and subjected to the loads 
a^ x (ai, s)n(aj, s) on S s . Therefore, by derivation with respect to s', we can 
prove that d/dj! {Cg" x ) is the elastic stress field in fi s , resulting from the 
imposed displacement on cMl and the loads d/d-f! (a < tJ lx r\) (oti,s) applied 

on S s . The equation l|43(l shows that these loads are equal to — h div <r m (ai, s). 

By our choice of the coordinate system we have the equality h — d r(s) /ds. 
For our purposes it is not restrictive to assume that when s tends to u s tends 
to u and e(u s ) tends to e(u). Indeed, there is a positive constant M such that 
for all s < 1 

c f I Vu s | 2 dss + W^HSu.) < / w{s(u s )) dx + W^^SuJ < M 
Jn Jn 

where c is the ellipticity constant of the tensor C. The last inequality is true 
because for all s the elastic energy of u s is smaller that the elastic energy of u 
and Sy; C K U dB r ^ s )(dK). By the compactness theorem of Bcllcttini, Coscia 
& Dal Maso [BCDM] and the regularity of u s , there exists a sequence Sh — > 
such that u~ Sh converges punctually to u and £(u~ Sh ) converges punctually to 
e(u). 

Therefore when s tends to a m (ai,s) = a s (ai,s) tends to a(ai,s), hence 
div a m tends to (because div a = Q). The conclusion is that — h div a m (ai, s) 
tends to too, so the displacement solution of the problem 

div cr(Vv) = in Cl s 

v = on cT2 

(j(Vv)n = on (both sides of) K n £l s 

cr(Vv)n = — h div <j m (ai,s) on S s , 
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which is d/dji / (u™), tends to 0. In particular 



lim H = _ JL u «« = o . (44) 

' ly=o 



We remark that 



o\ d w\vu s ax - U. W VU QI 1 /■ 

hm M ^ M MB - (a) = ilim / (on)-- 

S^O S 2 S — Jg 

This equality together with 1441) concludes the proof of (|38l) 



2 implication. We want to prove now that if 

u g ^ 2 '°°(u°) n w^(n \ k, R n ) 

and for any open set B C il 

I X2 I (u)(B) < CM+(u) c (B) , 

then u = u(if, u°). 

The regularity assumptions concerning u and K allow us to make the further 
calculations. As previously, we use the notations: a = cr(Vu), w — io(Vu). 

Let us consider a vector field r\ satisfying i|32fl. Then 

fc2(u,7?) = aijUi.k-rjk.j - w divq dx = / -gu^vli^ T]k - [w r)i - crijU ltk r]k] 4 dx . 

Jn 

We deduce that the absolute continuous part of the measure | K2 | (u) has the 
density | 07^11^ |. The hypothesis (|31[) and the decomposition l|29|l imply that 
this density is equal to 0, hence: 

div er(Vu) = a.e. in O . 

Therefore a is an admissible stress in the sense of the definition 3.1.1. . The 
before mentioned result of Del Piero [DP] implies that crn has no jumps over 
K. 

The measure | K2 | has a part concentrated on K, i.e. the absolute contin- 
uous part of I K2 | with respect to the measure H,™^ 1 . The density of this part 
can be calculated from the expression of fe2: 

fc2(u, 77) = - [w rji - aijUi^Vk] i d.T 
Jn 

Indeed, 1321) and the regularity assumptions over K and u allow us to use a 
flux-divergence formula in order to prove that the density of this part of | K2 \ 
is I crn |. Again, (13 If) and the decomposition (|29|l imply that this density is 
equal to 0, hence: 

crn = on (both sides of) K 
This concludes the proof. • 
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8 Final considerations 



Let us return to relation 1|36[1 from the proof of theorem 6.1., first implication. 
We have considered there u = u(K, u ) and r\ an arbitrary vector field, such 
that 

I 77 |< 1 (45) 
r\ ■ n = on K . 

This field generates the one-parameter flow </> s , which represents an arbitrary 
smooth crack propagation curve with sub-unitary initial velocity 77. The energy 
release rate, defined by l|17fl . associated to this curve, at the moment s = will 
be denoted by £ It has the expression: 

£(77) = - lim -^ w ( Vu ^ s ' u °)) dx Z -/n w ( Vu ) dx 

s—>0 s 

With this notation the relation ll.'UiH becomes: 

-£{fj) < -K2(l Q ,r),u°) . 
In fact we have proved more, namely that 

-CM + (uf(n) < -£( V ) < -K2(l n ,i h u°) . 
Therefore we have the following relation: 

I K2 I (u(K,u°))(n) < sup {£(77) : 77 satisfies (gSJ } < CM + (u(K, u°)) cr (0) . 

(46) 

This means that the supremum of the initial energy release rate, taken over all 
possible smooth propagations of the initial crack K with sub-unitary speed, is 
bounded by the concentration coefficients | K2 | and CM + , calculated for the 
displacement of the body with the initial crack K . 

The assumption that Griffith and Irwin criterions of brittle crack propaga- 
tion are equivalent is usual in classical fracture mechanics. From our considera- 
tions it is straightforward that this assumption means that at any moment the 
coefficients | K2 \ and CM + are equal, which leads to the following physical 
statement: during the evolution, the crack takes energy from the body with 
maximal speed. We leave for further examinations whether this statement can 
be generally proved or it has to be supposed. 



August 1, 1997 Marius Buliga 

Institute of Mathematics 
of the Romanian Academy 
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